Land cover dataset of the China Central-Asia West-Asia Economic Corridor from 1993 to 2018.

Land Cover (LC) maps offer vital knowledge for various studies, ranging from sustainable development to climate change. The China Central-Asia West-Asia Economic Corridor region, as a core component of the Belt and Road initiative program, has been experiencing some of the most severe LC change tragedies, such as the Aral Sea crisis and Lake Urmia shrinkage, in recent decades. Therefore, there is a high demand for producing a fine-resolution, spatially-explicit, and long-term LC dataset for this region. However, except China, such dataset for the rest of the region (Kyrgyzstan, Turkmenistan, Kazakhstan, Uzbekistan, Tajikistan, Turkey, and Iran) is currently lacking. Here, we constructed a historical set of six 30-m resolution LC maps between 1993 and 2018 at 5-year time intervals for the seven countries where nearly 200,000 Landsat scenes were classified into nine LC types within Google Earth Engine cloud computing platform. The generated LC maps displayed high accuracies. This publicly available dataset has the potential to be broadly applied in environmental policy and management.

Measuring and evaluating SDG indicators with Big Earth Data.

The United Nations 2030 Agenda for Sustainable Development provides an important framework for economic, social, and environmental action. A comprehensive indicator system to aid in the systematic implementation and monitoring of progress toward the Sustainable Development Goals (SDGs) is unfortunately limited in many countries due to lack of data. The availability of a growing amount of multi-source data and rapid advancements in big data methods and infrastructure provide unique opportunities to mitigate these data shortages and develop innovative methodologies for comparatively monitoring SDGs. Big Earth Data, a special class of big data with spatial attributes, holds tremendous potential to facilitate science, technology, and innovation toward implementing SDGs around the world. Several programs and initiatives in China have invested in Big Earth Data infrastructure and capabilities, and have successfully carried out case studies to demonstrate their utility in sustainability science. This paper presents implementations of Big Earth Data in evaluating SDG indicators, including the development of new algorithms, indicator expansion (for SDG 11.4.1) and indicator extension (for SDG 11.3.1), introduction of a biodiversity risk index as a more effective analysis method for SDG 15.5.1, and several new high-quality data products, such as global net ecosystem productivity, high-resolution global mountain green cover index, and endangered species richness. These innovations are used to present a comprehensive analysis of SDGs 2, 6, 11, 13, 14, and 15 from 2010 to 2020 in China utilizing Big Earth Data, concluding that all six SDGs are on schedule to be achieved by 2030.

Construction of a Territorial Space Classification System Based on Spatiotemporal Heterogeneity of Land Use and Its Superior Territorial Space Functions and Their Dynamic Coupling: Case Study on Qionglai City of Sichuan Province, China.

Territorial space classification (TSC) provides the basis for establishing systems of national territory spatial planning (NTSP) and supervising their implementation in China, thus has important theoretical and application significance. Most of the current TSC research is related to land use/land cover classification, ignoring the connection of the NTSP policies and systems, failing to consider the spatiotemporal heterogeneity of land use superior territorial space functions (TSFs) and the dynamic coupling between land use and its superior TSFs on the result of TSC. In this study, we integrated the factors influencing the connection of NTSP policies and systems and established a theoretical framework system of TSC from the perspective of spatial form and functional use. By integrating the q-statistic method with spatiotemporal geographical analysis, we propose a method to construct a TSC system for Qionglai City of Sichuan Province in China based on the spatiotemporal heterogeneity of land use superior TSFs and the dynamic coupling between land use and its superior TSFs. It makes up for the deficiency of directly taking land use/land cover classification as TSC and solves the problems of ignoring the spatiotemporal heterogeneity of land use superior TSFs and the dynamic coupling between land use and its superior TSFs. Using this method, we found that the TSC of Qionglai City consists of 3, 7, and 14 first-, second-, and third-level space types, respectively. The key findings from this study are that land use superior TSFs show spatiotemporal heterogeneity in Qionglai, and coupling effects in spatial distribution were noted between land use types and their superior TSFs, as was temporal heterogeneity in the coupling degree and the structure of the TSFs corresponding to the land use types, which show obvious dynamics and non-stationarity of the functional structure. These findings confirm the necessity of considering the spatiotemporal heterogeneity of land use superior TSFs and the dynamic coupling between land use and its superior TSFs in TSC. This method of establishing a TSC system can be used to address a number of NTSP and management issues, and three examples are provided here: (a) zoning of urban, agricultural, and ecological space; (b) use planning of production, living and ecological space; (c) delimitation of urban development boundary, permanent basic farmland protection redline, and ecological protection redline.

Assessment of five satellite-derived LAI datasets for GPP estimations through ecosystem models.

Ecosystem models have been widely used for obtaining gross primary productivity (GPP) estimations at multiple scales. Leaf area index (LAI) is a critical variable in these models for describing the vegetation canopy structure and predicting vegetation-atmosphere interactions. However, the uncertainties in LAI datasets and the effects of their representation on simulated GPP remain unclear, especially over complex terrain. Here, five most popular datasets, namely the Long-term Global Mapping (GLOBMAP) LAI, Global LAnd Surface Satellite (GLASS) LAI, Geoland2 version 1 (GEOV1) LAI, Global Inventory Monitoring and Modeling System (GIMMS) LAI, and Moderate Resolution Imaging Spectroradiometer (MODIS) LAI, were selected to examine the influences of LAI representation on GPP estimations at 95 eddy covariance (EC) sites. The GPP estimations from the Boreal Ecosystem Productivity Simulator (BEPS) model and the Eddy Covariance Light Use Efficiency (EC-LUE) model were evaluated against EC GPP to assess the performances of LAI datasets. Results showed that MODIS LAI had stronger linear correlations with GLASS and GEOV1 than GIMMS and GLOMAP at the study sites. The GPP estimations from GLASS LAI had a better agreement with EC GPP than those from other four LAI datasets at forest sites, while the GPP estimations from GEOVI LAI matched best with EC GPP at grass sites. Additionally, the GPP estimations from GLASS and GEOVI LAI presented better performances than the other three LAI datasets at crop sites. Besides, the results also showed that complex terrain had larger discrepancies of LAI and GPP estimations, and flat terrain presented better performances of LAI datasets in GPP estimations. Moreover, the simulated GPP from BEPS was more sensitive to LAI than those from EC - LUE, suggesting that LAI datasets can also lead to different uncertainties in GPP estimations from different model structures. Our study highlights that the satellite-derived LAI datasets can cause uncertainties in GPP estimations through ecosystem models.

[An Improved DDV Method to Retrieve AOT for HJ CCD Image in Typical Mountainous Areas].

Domestic HJ CCD imaging applications in environment and disaster monitoring and prediction has great potential. But, HJ CCD image lack of Mid-Nir band can not directly retrieve Aerosol Optical Thickness (AOT) by the traditional Dark Dense Vegetation (DDV) method, and the mountain AOT changes in space-time dramatically affected by the mountain environment, which reduces the accuracy of atmospheric correction. Based on wide distribution of mountainous dark dense forest, the red band histogram threshold method was introduced to identify the mountainous DDV pixels. Subsequently, the AOT of DDV pixels were retrieved by lookup table constructed by 6S radiative transfer model with assumption of constant ratio between surface reflectance in red and blue bands, and then were interpolated to whole image. MODIS aerosol product and the retrieved AOT by the proposed algorithm had very good consistency in spatial distribution, and HJ CCD image was more suitable for the remote sensing monitoring of aerosol in mountain areas, which had higher spatial resolution. Their fitting curve of scatterplot was y = 0.828 6x-0.01 and R2 was 0.984 3 respectively. Which indicate the improved DDV method can effectively retrieve AOT, and its precision can satisfy the atmospheric correction and terrain radiation correction for Hj CCD image in mountainous areas. The improvement of traditional DDV method can effectively solve the insufficient information problem of the HJ CCD image which have only visible light and near infrared band, when solving radiative transfer equation. Meanwhile, the improved method fully considered the influence of mountainous terrain environment. It lays a solid foundation for the HJ CCD image atmospheric correction in the mountainous areas, and offers the possibility for its automated processing. In addition, the red band histogram threshold method was better than NDVI method to identify mountain DDV pixels. And, the lookup table and ratio between surface reflectance between red and blue bands were the important influence factor for AOT retrieval. These will be the important research directions to further improve algorithm and improve the retrieve accuracy.